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Abstract 

The cross section for the diffractive reaction y* + p ^ cc + p with 
a real or virtual photon is calculated in the nonperturbative two-gluon 
exchange model of Landshoff and Nachtmann. Numerical predictions 
are given for cross sections and spectra at typical HERA values of c.m. 
energy and photon virtuality. The contribution of charm to the diffrac¬ 
tive structure function is evaluated and found to be rather small in the 
model, and the ratio between the production rates for bb and cc is tiny. 
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1 Introduction 


The discovery of diffractive events in ep collisions at HERA [jl[ has triggered a 
large amount of experimental and theoretical work and greatly increased our 
knowledge of the physics of diffraction, or of the pomeron. Several models 
give a reasonable description of the data at the time being, but we are far yet 
from a coherent picture of the mechanisms at work in terms of QCD. One can 
hope that the detailed study of the diffractive hnal state will lead to further 
progress in this direction. Charm production looks promising in this respect, as 
predictions for this process differ widely between various models ii i ii- 
In this paper we use the approach due to Landshoff and Nachtmann (LN) 
to model the soft pomeron by the exchange of two nonperturbative gluons. 
We present differential cross sections for the diffractive dissociation of a real or 
virtual photon into a cc-pair. Provided that the invariant mass of the diffractive 
hnal state is not too large its cc-component should give a fair approximation of 
inclusive diffractive charm production. In the following section we give some 
details of the model and of the calculation, in sec. ^we present our results, and 
in sec. ^ we summarise our hndings. 

2 Diffractive cc-production in the LN model 

The Landshoff-Nachtmann model has been introduced and described in ||^, 
It approximates the soft pomeron by the exchange of two nonperturbative 
gluons, with a propagator —g^^D{p) instead of the perturbative in 

Feynman gauge. In several processes one can express the scattering amplitude 
in terms of a few moments of the function D{1‘^) and thus does not need to 
know its detailed form. Here we only need the integral 

= , 1 ) 

where /do ~ 2.0GeV~^ and po ~ 1.1 GeV have been extracted from experi¬ 
mental data p, H. The parameter pg provides the characteristic scale for the 
dependence of D{p) on the gluon virtuality /^, and stands for the strong 
coupling in the nonperturbative region which dominates the /^-integration in 
(]^). It will be taken as 1 here 0. 

We now apply this model to the reaction 

7 * + p ^ cc -f p . (2) 

In the following we use the conventional variables x,y, s,t,Q‘^,W‘^ for deep 
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Figure 1: Two of the four Feynman diagrams contributing to the imaginary 
part of the amplitude for p + 7 * p + cc. The other two are obtained by 
reversing the charge flow of the upper quark line. The lower line stands for 
a constituent quark in the proton as explained in O 
denote nonperturbative gluons. 


and the dashed lines 


inelastic scattering, M for the invariant mass of the cc-pair, and 


/3 




+ m2 -1 

hF2 + Q2 


(3) 


We denote with Pt the transverse momentum of the charm quark with respect 
to the photon momentum in the 7 *p c.m. 

In the high-energy limit the scattering amplitude for our process is domi¬ 
nated by its imaginary part and we can use the cutting rules to calculate it. 
Then the diagrams contributing to (^) are those in £g. and the ones obtained 
by reversing the charge flow of the upper quark line. In each diagram there 
is one off-shell quark, the characteristic scale for its virtuality being given by 


2 Pt + < 
1-/3 


(4) 


For charm production the large quark mass me protects this quark from be¬ 
coming infrared so that a perturbative treatment of the quark sector should be 
safe, even in the photoproduction limit. 

The cross section for photons with transverse or longitudinal polarisation 
in the 7 *p frame has been calculated in [O, 


T^. It reads 
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where Cc = 2/3 is the electric charge of the charm quark in units of the positron 
charge, and Fi{t) the Dirac form factor of the proton. We have approximated 
f = 0 for the squared momentum transfer from the proton, except in Fi{t) 
and a]p(t), which accounts for most of the f-dependence in the cross section 
|T2| . apit) ^ 1.085+ f/(2 GeV)^ is the soft pomeron trajectory as observed in 
hadronic reactions [|T3. It is introduced by hand in the LN model in order to 


make contact with experiment; the approximation of bare two-glnon exchange 
(fig. i would give a factor instead of ^pd)) the cross section (^. 
We thus assume that the energy dependence of diffractive charm production is 
given by the soft pomeron, and fnrthermore that ^ is the correct dimensionless 
variable to be raised to the Regge power (l/,^)“^h) the amplitnde. In our 
numerical applications we will impose an upper cut of .^ < 0.05 to remain 
in a region where the exchange of a pomeron dominates that of secondary 
trajectories. 

The expressions 


= 1-2 
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in (ID contain integrals Li{P^,w) and L 2 {PY,w) over the virtnality of the 
exchanged glnons. 


roo ^ ^ 

Li{Ppw)= / dl‘^[af^D{-iyf fi{v,w) , 
Jo 
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with 
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f2{v,w) = 1 
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( 9 ) 


A simple way to approximate Li is to Taylor expand the function fi{v,w) 
about n = 0 and to keep only the leading term 


fi{v,w) V ■ 


dfi{v,w) 


dv 


v=0 


( 10 ) 


after which the integrals reduce to the moment (|l|) of the gluon propagator. 
This is however not very good for small where the variable v becomes 
V = {1 — P) ■ P /nil needs a good approximation of the functions 

fi{v), f 2 {v) for V from zero to order one. A better approximation which also 
leads to the moment i) is 


fi{v,w) ^ 


V 

Vo 


fi{vo,w) , 


.0 = ^ 


( 11 ) 


with /q ~ Pq. Whereas in (pTl) one approximates the curve fpv) by its tangent 
at n = 0 the approximation (^) uses instead the line that intersects the curve 
at n = 0 and v = vo, the corresponding range in P from 0 to /g being the 
dominant region of integration in Lj. We have varied the parameter /g between 
/Uq/2 and 2/iQ and found variations of up to a factor 1.6 in the integrated 
cross sections ut and and variations of less than a factor 2.1 for cross 
sections differential in or in together with some change in the shape of 
the spectra. The effects are stronger at small or zero and more pronounced 
in the P^ - than in the M^-spectra. These variations may be seen as reflecting 
our uncertainty about the exact shape of the nonperturbative gluon propagator, 
which determines the value of /g for which the approximation ( pd^ is best. As a 
benchmark we have compared the approximated integrals with the exact ones 
for the model gluon propagator used in [^, 


D{-P) oc 


1 + 


P 


(n - l)fil 


n > A , 


( 12 ) 


where the proportionality constant can easily be obtained from (|l|). With ([Tl|) 
the errors of the approximation stay below 10% for all P^ and w we need, 
whereas the tangent approximation (0) has errors of 50% and more when one 
goes to P| = 0. 

Finally a comment is in order about the value of the strong coupling in our 
calculation. In the integral (|^) it is taken at a nonperturbative scale given by 

^The definitions of v and w here differ from those in im. 
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the dominant virtuality of the exchanged gluons. However, the hrst gluon in 
all diagrams couples at its upper end to an off-shell quark whose virtuality is 
in the perturbative region. We choose to take the coupling for this vertex at 
the scale in all four diagrams, and at the small gluonic scale for the other 
three vertices. One might argue that is the typical scale for the entire upper 
parts of the diagrams, and that one should also take the coupling of the second 
gluon to the upper quark line at this scale. Note however that, since we use 
the cutting rules, the parts of the diagrams to the left and the right of the cut 
lines can be considered independently, and that the part to the right of the 
cut is just on-shell quark-quark or quark-antiquark scattering with no large 
virtuality involved. We are aware though that our choice is only a guess. 

It is clear that this question of scales leads to an uncertainty in the nor¬ 
malisation of our predictions. The main problem is not so much whether A^ is 
the best choice of a hard scale, which is a common problem of all leading order 
calculations. Since the moment (|^) includes the nonperturbative coupling 
we have to multiply the cross section with as(A^)/a^°^ if at one of the four 
quark-gluon vertices we take the perturbative coupling, is not well con¬ 
strained by phenomenology or theory and other choices than = 1 which 
we adopt here have indeed been made [|^ . 


3 Results 


We will now give some numerical predictions for the diffractive production 
of charm. We stress once again that we calculate the cross section for the 
diffractive final state being a quark-antiquark pair, which does not include 
events with cc and additional gluons at parton level. It is also different from 
cc-production through photon-gluon fusion where the gluon is a parton emitted 
by the pomeron in the description of Ingelman and Schlein and where 
the hnal state contains a pomeron remnant. Finally it excludes events in 
photoproduction where the photon is resolved. 

Let us start with photoproduction. For W = 220 GeV and ^ < 0.05 we 
obtain a total rate cr( 7 p —>■ ccp) = 57 nb. Spectra in and are shown in 
fig- 1 and The -spectrum is well described by a power behaviour 


da. 


IP 


dP^ 


oc (Hp + mJ 


2\-5 


(13) 


with an exponent 6 between 3.7 and 4.3 in the -range of fig. The spectrum 
of diffractive mass behaves approximately like da.yp/dM‘^ cx to the right 

of its peak. 
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Figure 2: Spectra in and in for 7 *p —ccp at PF = 220 GeV with a 
cut ^ < 0.05. The curves are for photoproduction and for electroproduction at 
= 20 GeV^ with transverse or longitudinal photons. 
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Figure 3: The spectra of £g. § (a) for a wider range in P^, as a function of 
P^ + rril in a double logarithmic plot. For photoproduction one recognises the 
power behaviour (|T^. 
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Let us move on to electroproduction and first focus on the dependence on 
the 7 *p c.m. energy W of the ep cross section, obtained from the usual relation 

dcTep = ^ ^ {(1 - 2/ + 2 /V 2 ) dar + (1 - 2/) dax} (14) 

where we have used the approximation x -C 1. It is determined by different 
effects: 

1. the integration element dW"^/W"^ in (0) 

2. the ^/-dependent factors multiplying dax and dax, note that at x -C 
1 one has y ~ W'^/s. They decrease with W. As they are different 
for transverse and longitudinal photons one cannot calculate their effect 
without knowing the relative contribution of ax and ax- This problem is 
of course absent in the photoproduction limit. 

3. the cut in ^ which ensures that the selected events are dominated by 
pomeron exchange. It leads to an upper limit on the diffractive mass M 
that depends on W and and whose importance is strongest for small 
W and large Q^. This effect can be circumvented by using a hxed cut on 
M chosen such that ^ is always small enough, but at the expense of the 
total rate used in the analysis. 

4. the dependence on W of dax^x/{dP^dM'^ dt). This gives direct informa¬ 

tion about whether or not this process is dominated by the soft pomeron. 
In our model it comes from the factor ^^(i-ojpp)) with -C 

it is which is quite flat given that the soft pomeron intercept 

is close to 1. 

We find an ep cross section of 120 pb for y/s = 296 GeV, ^ < 0.05, integrated 
over Q'^ from 7.5 GeV^ to 80 GeV^ and W from 50 GeV to 220 GeV. Table 
(a) gives the ep cross section in three different bins of W. They are spaced 
logarithmically to take out the trivial effect of point 1. above. One might 
also choose the binning to include the factor 1 — y + |/^/2 if one assumes the 
longitudinal contribution to the cross section to be small; one then directly 
extracts the fV-dependence of the integrated 'y*p cross section ax- The effects 
of points 2. and 3. are responsible for the decrease of the cross section from 
the second to the third IV-bin in our numerical example. Apart from this one 
sees however clearly the flat behaviour in W characteristic of soft pomeron 
exchange. It might be useful to analyse experimental data in this way: it 
focuses on the fV-dependence of the 7 *p cross section and apart from points 
2. and 3. does not involve the details of its dependence on or Q"^, so it 
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is relatively model independent. Also it requires only binning in one variable 
and thus makes best use of the available statistics which is likely not to be 
abundant at HERA. If there were a strong departure from the soft pomeron 
energy dependence for these events it should be seen in such an analysis. 


Table 1: (Jep for ep —>• epcc with -y/s = 296 GeV and ^ < 0.05. (a) For from 
7.5 GeV^ to 80 GeV^ and three logarithmically spaced bins in W. (b) For W 
from 50 GeV to 220 GeV and three logarithmic bins in Q^. 


(a) 



(6) 



50 to 82 

1U[ GeV] 
82 to 134 

134 to 220 

7.5 to 16.5 

Q2[ GeV^j 
16.5 to 36.3 

36.3 to 80 

39 pb 

44 pb 

40 pb 

73 pb 

36 pb 

14 pb 


In table (6) we give the ep cross section in logarithmically spaced Q^-bins. 
Since y is nearly independent of at small x this directly shows the Q^- 
dependence of the weighted sum of the integrated transverse and longitudinal 
cross sections. 

Examples of ep spectra are shown in fig. ^ and 7 *p spectra for transverse 
and longitudinal photons in fig. In the M^-spectrum of fig. |] (6) we can 
see that the contribution from longitudinal photons is much smaller than from 
transverse ones except at very small M^. This is explained by the factor (77^ + 
m^)/M‘^ in the differential cross section for longitudinal polarisation, cf. eq. (^. 
The -spectrum for transverse photons is less steep than in photoproduction 
at its lower end, but at large it becomes similar in slope and normalisation. 
The dip in the P^ -spectrum for longitudinal polarisation is a consequence of a 
zero in the longitudinal cross section at fixed P^ and which is due to the 
behaviour of the integral L 2 {Px, w) in ([^. Unfortunately it becomes completely 
swamped in the ep spectrum by the transverse contribution, cf. fig. (a), so 
that this effect is unlikely to be observed. 

The M^-spectra can be rewritten in terms of the charm contribution to the 
diffractive structure functions and dehned by 


da, 


ep 


Airat 


dx dQ'^ d^ dt 


xQ‘^ 


(l-y + yy2)F^^^^ + {l-y)F[ 


D(4) 


Fr 


D(4) 


= Fr 


D(4) 


+ Ff 


(15) 


where T and L stand for the contributions of transverse and longitudinal pho¬ 
ton polarisation as usual. These functions depend on the kinematic variables 
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Figure 4: Spectra in and in for ep —*• epcc at a/s = 296 GeV, integrated 

over ^ < 0.05, PF = 50 GeV to 220 GeV and = 7.5 GeV^ to 80 GeV/ 


or, equivalently, on In models with Regge factorisation 

this dependence factorises into 

(?. /3. Q". *) = fA(, t) ■ Ffpn, Q\ t) , (16) 

where in the partonic interpretation of Ingelman and Schlein f]p{^, t) gives 
the pomeron flux from the proton and Q^, t) are the structure functions 

of the pomeron. 

Our model has this factorisation property, but this is because it is put in 
rather than being one of its predictions: We calculate two-gluon exchange to 
leading order in so that [3, t) depends on ^ via a global factor, 

and then we modify the exponent of ^ by hand introducing the soft pomeron 
trajectory, which preserves factorisation. Using the flux factor 

= ( 17 ) 

we obtain the curves shown in fig. ^ for the charm structure functions (cc) 
and Ff/’{cc) of the pomeron. We observe that the transverse structure function 
shows Bjorken scaling for rather small fd. It can be shown from the expressions 
d^) and (H) that F^ (cc) scales under the condition that the integrated cross 
section is dominated by values of P^ that are small compared with the upper 
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Figure 5: Charm structure functions (cc), F[’{cc) of the pomeron for trans¬ 
verse and longitudinal photon polarisation at different values of Q^. 


kinematical limit M^/4 — m^. Because with our values of the diffractive 
mass M is not far from the production threshold 2mc for moderate and large 
(3 this condition cannot be satished in this region and there is no scaling. Note 
also that at the value of j3 corresponding to M = 2mc both F^(cc) and F^(cc) 
vanish. 

We wish to compare the rate of charm production with the inclusive diffrac¬ 
tive cross section for light flavours u, d, s, which can also be calculated in our 
model O, ro. Let us remark that for light flavours the quark virtuality 


in dJ) can become small so that one has to assume that a perturbative treat¬ 
ment of the quarks, possibly with an effective quark mass of some 100 MeV, 
is good enough for this process. We hnd that actually the best description of 
the HERA data with our model is obtained when just taking current masses 
for the M, d and s quarks , which we have done in the results to be shown. 
For light quarks the approximations (TO) or (P!TD of the loop integrals Li, L 2 
can no longer be used unless is large so that one has to resort to a specihc 
form of the gluon propagator D(P). The results we show have been obtained 
using the model propagator (0) with n = 4 and freezing the running coupling 
as(A^) in (^) when it becomes equal to 1. Frf shows scaling over the entire 
/9-range since there is no strong threshold effect as for charm at large f3. The 
longitudinal structure function F[’ is found to behave roughly like 1/Q^ at 
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Figure 6: HERA data for the ^integrated diffractive structure function 
Q^) — /^Tir /3, g^, t) compared with the result in the LN 

model. Also shown is the longitudinal contribution Fj^ to F^. Integration 
limits are ^min = 6 . 3 - 10 “"^, ^max = 0.01 in the case of ZEUS and (min = 3 - 10 “"^, 
(max = 0.05 for the HI data. 


hxed j3 and only gives a signihcant contribution to F^ when (3 is large. 

Fig. I shows the result of our calculation of the ^integrated diffractive 
structure function F^{l3,Q^) = d( J dt F^^'^\(, l3,Q‘^,t) together with 
HERA data m We only show one per experiment as the dependence 
of F^^^\(, (3^ g^, t) on g^ is found to be weak in the data, in good agreement 
with our results. Remember that our calculation is at Born level and does not 
incorporate the effects of QGD evolution of structure functions. 

Regarding the overall normalisation we hnd that agreement is not too bad 
given that we have done a leading order calculation and taking into account the 
uncertainty in its normalisation due to the strong coupling at different scales 
we discussed in sec. 0. We also stress that the parameters of our model, /3g, /Xq 
and have all be determined from pre-HERA data and that in this sense 
our prediction is parameter free. As to the shape in (3 the data clearly do not 
show a decrease at small (3 as our result does. This is not surprising since 
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Figure 7: Pomeron structure function F^{cc) for charm compared with 0.05 
times for the three light flavours, (a) for = 16 GeV^ and {h) for 
g2 = 25 GeV^ 


we only calculate the gg-component of the diffractive dual state, and at small 
values of f3, i.e. at large diffractive mass M final states with additional gluons 
are expected to be dominant. Fig. ^ indicates that this might be the case for 
values of (3 well above 0.1. We remark that the more recent ZEUS data ITS 


indicate a rise of as {3 becomes small whereas the preliminary HI data 

give a very flat behaviour over the entire /d-range. 


In fig. 1^ we compare the predictions of the model for FI 


Fi 


cc 


and F^ for the 


three light flavours, keeping in mind that neither is expected to be a complete 
description at small (3. The curves for light quarks are scaled down by a factor 
of 20. The fraction of charm comes out quite small,^ it is not larger than 
5% and decreases with (3. For (3 > 0.4, where we expect the gg-component to 
dominate the final state, we find the fraction of charm in F^ to be below 4%. 

We observe that for the longitudinal pomeron structure function the con¬ 
tribution of charm can be considerable. It is largest for (3 around 0.5 where 
the ratio between FF(cc) and 


F^ for all flavours can reach values close to 


^The results for the charm contribution to the diffractive structure function presented 
here are significantly smaller than those given in because there we used the running 
coupling as{Px) which unlike as(A^) is not limited from above by as(rn^) and was frozen 
when it reached the value 1. 
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0.6. This is larger than the ratio of the squared electric charges so that apart 
from the charge the cross section for charm exceeds that of a light gg-pair. 
Note that the differential longitudinal cross section has a suppression factor 
{P^ + m^)/M‘^ at small which is less efficient for large quark mass This 
phenomenon should however be difficult to observe since the longitudinal part 
of is tiny at f3 0.5 and only visible at rather high f3 and not too large 
where is below the charm threshold. 

To conclude we wish to report an observation regarding the ratio of the rates 
for diffractive cc- and &&-production. We compare the integrated cross sections 
for photoproduction, so that the scale cannot influence this ratio, taking 
W = 220 GeV and imposing ^ < 0.05. We checked that at this W the effect 
of the phase space reduction through the cut in ^ is about the same for charm 
and bottom: calculating the cross sections for very high W where the .^-cut has 
no effect we verified that the increase of the cross sections is almost entirely 
due to the factor from pomeron exchange. With nic = 1.5 GeV and 

mh = 4.5 GeV we find that the cross section for bottom production is about 
440 times smaller than for charm, implying that this process should be quite 
impossible to observe at HERA. Taking out the squared electric charges and 
the running strong coupling, which for simplicity we take here at fixed scales 
ml or ml, we obtain 

a( 7 P ^ ccp) el as{ml) _ 

o'i'jp^bbp) elasiml) \mcJ 

Apart from the effect of the running coupling the integrated cross section ap¬ 
pears to scale with the quark mass approximately like l/rrig. To check this we 
have calculated the photoproduction cross section as a function of mq between 
mq = 1.5 GeV and 15 GeV, for very large W so that the .^-cut has no effect. 
Dividing by the running coupling as{ml) we find indeed an approximate power 
behaviour in 1/mq with an exponent between 3.6 and 4.2. A behaviour in l/m^ 
looks like the effect of the off-shell propagators in the amplitude whose denom¬ 
inators are limited by m^. Notice that the numerators of the propagators, 
which also contain one power of mq do not seem to enter in the same way, their 
role is more complicated because the numerator of a quark propagator has a 
Dirac matrix structure. 


4 Summary 

We have calculated diffractive production of a cc-pair in 7 *p collisions with real 
or virtual photons in the model of nonperturbative two-gluon exchange due to 
Landshoff and Nachtmann. This allowed us to give numerical predictions for 
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diffractive charm production at HERA, for cross sections and spectra in 
and P^. 

In photoproduction we hnd a 'jp cross section in the region of 60 nb for W 
around 200 GeV. The mass spectrum peaks at rather low values of M and then 
falls off roughly like 1/M^, whereas the spectrum of the transverse momentum 
approximately behaves like a power of with an exponent around —4. 

In diffractive DIS the ep cross section we obtain is of order 100 pb for 
from 7.5 GeV^ to 80 GeV^ and W from 50 GeV to 220 GeV. The 7 *p cross 
sections exhibit a flat dependence on W, typical of models with soft pomeron 
exchange. We suggest that even a coarse logarithmic binning in W of the ep 
cross section should be useful to test this in the data. 

Expressing the cc mass spectra in terms of the diffractive charm structure 
functions we hnd that the transverse contribution F^^‘^\cc) does not scale for 
f3 > 0.3 due to the restricted phase space at HERA values of Q^. The longi¬ 
tudinal structure function F^^^\cc) for charm is comparable to the transverse 
one at small diffractive mass M, i.e. close to the largest kinematically allowed 
f3 at given at lower f3 it is negligible. 

We have then compared F 2 ^^\cPj with the diffractive structure function 
F 2 ^^'^ for light havours, calculated in the same model, which reproduces the 
HERA data within a factor of 2 or so, except in the region of small (3 where 
the neglect of hnal states other that qq becomes a bad approximation. The 
contribution of cc to the diffractive structure function comes out below 

5% in our model. 

The strong quark mass dependence of the integrated qq cross section can 
be understood in a simple way from the denominator of the propagator for the 
off-shell quark in the Eeynman diagrams, which appears squared in the cross 
section. Its typical value is given by (A^) = ((7^) +m^)/{l — P), where (P^) is 
some average P^. For photoproduction of heavy quarks we hnd indeed a quark 
mass dependence of the cross section approximately like l/rUg, which together 
with the effect of the running strong coupling and the different quark charges 
leads to a ratio of bb- to cc-production of 1/440. 
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